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Numerical calculations have been performed to elucidate unconventional electronic transport prop- 
erties in disordered nanographene ribbons with zigzag edges (zigzag ribbons). The energy band 
structure of zigzag ribbons has two valleys that are well separated in momentum space, related to 
the two Dirac points of the graphene spectrum. The partial flat bands due to edge states make 
the imbalance between left- and right-going modes in each vaUey, i. e. appearance of a single chiral 
mode. This feature gives rise to a perfectly conducting channel in the disordered system, i.e. the 
average of conductance {g) converges exponentially to 1 conductance quantum per spin with in- 
creasing system length, provided impurity scattering does not connect the two valleys, as is the case 
for long-range impurity potentials. Ribbons with short-range impurity potentials, however, through 
inter-valley scattering, display ordinary localization behavior. Symmetry considerations lead to the 
classification of disordered zigzag ribbons into the unitary class for long-range impurities, and the 
orthogonal class for short-range impurities. The electronic states of graphene nanoribbons with 
general edge structures are also discussed, and it is demonstrated that chiral channels due to the 
edge states are realized even in more general edge structures except for armchair edges. 

PACS numbers: 72.10.-d,72.15.Rn,73.20.At,73.20.Fz,73.23.-b 



I. INTRODUCTION 

Graphene being for many decades only a domain to 
theoretical studies, has recently been fabricated with in- 
genious methods and has initiated intensive and diverse 
research on this system^ The honeycomb crystal struc- 
ture of single layer graphene consists of two inequivalent 
sublattices and results in a unique band structure for 
the itinerant 7r-electrons near the Fermi energy which 
behave as massless Dirac fermion. The valence and 
conduction bands touch conically at two nonequivalent 
Dirac points, called Kj^ and K— point, which form a 
time-reversed pair, i.e. opposite chiralityi^ The chirality 
and a Berry phase of tt at the two Dirac points pro- 
vide an environment for highly unconventional and fasci- 
nating two-dimensional electronic properties, such as the 
half-integer quantum Hall effect j^i^ the absence of back- 
ward scatteringf^ 7r-phase shift of the Shubnikov-de Haas 
oscillations 

The successive miniaturization of the graphene elec- 
tronic devices inevitably demands the clarification of 
edge effects on the electronic structures and electronic 
transport properties of nanometer-sized graphenes. The 
presence of edges in graphene has strong implications for 
the low-energy spectrum of the 7r-electrons.^'^ There are 
two basic shapes of edges, armchair and zigzag which 
determine the properties of graphene ribbons. It was 
shown that ribbons with zigzag edges (zigzag ribbon) 
possess localized edge states with energies close to the 
Fermi leveli^'^'^ These edge states correspond to the non- 



bonding configurations as can be seen by examining the 
analytic solution for semi-infinite graphite with a zigzag 
edge for which the wave functions of the edge states re- 
side on one sublattice only.^ In contrast, edge states are 
completely absent for ribbons with armchair edges. Re- 
cent experiments support the evidence of edge localized 
states. ^"'^^ Also, graphene nanoribbons can experimen- 
tally be produced by using lithography techniques'^ 

The electronic transport through zigzag ribbons shows 
a number of intriguing phenomena such as zero- 
conductance Fano resonances ji^ii^ vacancy configuration 
dependent transport, valley filtering, half-metallic 
conduction^^ and spin Hall effect.^'* It is also expected 
that the edge states play an important role for the mag- 
netic properties in nanometer-sized graphite systems, be- 
cause of their relatively large contribution to the density 
of states at the Fermi energyjii2ii^iSai2ii^i2^ 

Since the graphene nanoribbons can be viewed as a 
new class of the quantum wires, one might expect that 
random impurities inevitably cause the Anderson local- 
ization, i.e. conductance decays exponentially with in- 
creasing system length L and eventually vanishes in the 
limit of L — !■ oo. However, it was shown that carbon 
nanotubes with long-ranged impurities possess a per- 
fectly conducting channeli^ Also, recently, the present 
authors reported that nanographene ribbon with zigzag 
edges possess one perfectly conducting channel if the im- 
purity potentials are long-ranged, induced by electronic 
states which originate from the edge statesi^ Recent 
studies show that perfectly conducting channels can be 
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FIG. 1: (a) Structure of graphene zigzag ribbon. The disor- 
dered region with randomly distributed impurities lies in the 
shaded region and has the length L. The lattice constant is a 
and the ribbon width A*' is defined as the number of the zigzag 
chains. Randomly distributed blue circles schematically rep- 
resent the long ranged impurities, (b) Energy dispersion of 
zigzag ribbon with A'^ = 10. The valleys in the energy disper- 
sion near k = 27r/3a (fc = — 27r/3a) originate from the Dirac 
Kj^[K-)-]yoint of graphene. The red-filled (blue-unfilled) cir- 
cles denote the right (left)-moving open channel at the energy 
-Eo (dashed horizontal line). In the left (right) valley, the de- 
generacy between right and left moving channels is missing 
due to one excess right (left) -going mode. The time-reversal 
symmetry under the intra- valley scattering is also broken. 



stabilized in two standard universality classes. One is 
the symplectic universality class with an odd number 
of conducting channels r^i^^i^ and the other is the uni- 
tary universality class with the imbalance between the 
numbers of conducting channels in two propagating di- 
rectionsi^i^'2^ The symplectic class consists of systems 
having time-reversal symmetry without spin-rotation in- 
variance, while the unitary class is characterized by the 
absence of time-reversal symmetry.— 

In this paper, we study the disorder effects on the elec- 
tronic transport properties of graphene zigzag ribbons. 
The edge states play an important role here, since they 
appear as special modes with partially flat bands and 
lead under certain conditions to chiral modes separately 
in the two valleys. There is one such mode of oppo- 
site orientation in each of the two valleys of propagating 
modes, which are well separated in fc-space. The key re- 
sult of this study is that for disorder without inter-valley 
scattering a single perfectly conducting channel emerges 
introduced by the presence of these chiral modes. This 
effect disappears as soon as inter- valley scattering is pos- 
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FIG. 2: Schematic figure of the scattering geometry at K- 
point in zigzag ribbons, where a single excess right-going 
channel exists. Here ric = 0, 1, 2, • ■ • . On the contrary, the 
has a single excess left-going mode. 
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FIG. 3: L-dependence of (a) the average, {g), (b) the cor- 
responding variance, var{g), and (c) corresponding log-plot 
\n{g — 1) of dimensionless conductance for zigzag ribbon with 
— 10, d/a = 1.5 (no inter-valley scattering), uq — 1.0, and 
riimp. ~ 0.1. More than 9000 samples with different impurity 
configuration are included in the ensemble average. 



sible. Therefore the behavior is determined by the range 
of the impurity potentials. As a function of the impu- 
rity potential range a crossover from the orthogonal to 
the unitary universality class occurs which is connected 
with the presence or absence of time reversal symmetry 
(TRS). 

We organize this paper as follows. In Sec. 2, we briefly 
introduce the electronic states of graphene nanoribbons 
with zigzag edges. In Sec. 3, some peculiar features in the 
scattering matrix formulation on the graphene nanorib- 
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bons are explained, where the non-square form of the 
reflection matrix and the implications of the perfectly 
conducting channel in nanoribbons are presented. Also, 
the results of numerical calculation are shown. In Sec. 4, 
we discuss the electronic states of graphene nanoribbons 
with general edge structures, and demonstrate that chi- 
ral channels due to the edge states is realized even in the 
general edge structures except for the armchair edge. In 
Sec. 5, the universality class of the nanographene ribbons 
is discussed. The conclusion is presented in Sec. 6. 



II. ELECTRONIC STATES 

We describe the electronic states of nanographites by 
the tight-binding model 



(1) 



where 7i.j = —1 if i and j are nearest neighbors, and 
otherwise. \i) represents the state of the p^-orbital on site 
i neglecting the spin degrees of freedom. In the follow- 
ing we will also apply magnetic fields perpendicular to 
the graphite plane which are incorporated via the Peierls 
phase: 
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where A is the vector potential. The second term in Eq. 
([T|) represents the impurity potential, Vi — V{ri) is the 
impurity potential at a position r,j . 

In Figllja), the graphite ribbon with zigzag edges 
(zigzag ribbons ) is shown. We assume edge sites are 
terminated by H-atoms. The ribbon width N is defined 
by the number of zigzag lines. The length of disordered 
region is defined as L. Fig. \Vlh) depicts the energy band 
structure of zigzag ribbon for N — 10. The zigzag rib- 
bons are metallic for arbitrary ribbon width. The most 
remarkable feature is the presence of a partly flat band at 
the Fermi level, where the electrons are strongly localized 
near the zigzag edge. Each edge state has a non- vanishing 
amplitude only on one of the two sublattices, having, 
thus, non-bonding character. However, in a zigzag rib- 
bon of finite width, two edge states coming from both 
sides, have finite overlap. Because they are located on 
different sublattices, they mix into a bonding and anti- 
bonding configuration. In this way the partly flat bands 
acquire a dispersion and become conductive except at 
exactly E = 0. Note that the overlap is increasing as 
k deviates from ±7r/a, because the penetration depth of 
the edge states increases and diverges at fc = ±27r/3a, 
where a is the lattice constant. 

We briefly discuss here the relation between valleys in 
the zigzag ribbons and two-dimensional graphene. The 
electronic states near the Dirac point can be described 
by the k • p Hamiltonian 




FIG. 4: Impurity range, d/a, 
dimensionless conductance, (g) 



dependence of the averaged 
for zigzag ribbon with N — 



10, uo = 1.0, L/a = 5000 and rn, 



0.1. More than 1000 



samples with different impurity configuration are included in 
the ensemble average. 



acting on the 4-component pseudo-spinor Bloch functions 
^{r) = [(f>K+A,(t>K+B,'pK_A,(l)K_B], which character- 
ize the wave functions on the two crystalline sublattices 
(A and B) for the two Dirac points (valleys) K±. Here, 7 
is the band parameter, kx(ky) are wavenumber operators, 
and r° is the 2x2 identity matrix. Pauli matrices cr^'i''^ 
act on the sublattice space [A, B), while r^'^'^ on the 
valley space {K±). Since the outermost sites along 1"* 
(A^*'*) zigzag chain are B(A)-sublattice, an imbalance be- 
tween two sublattices occurs at the zigzag edges leading 
to the boundary conditions 



(4) 



Hk.p = l hicy"" ® T°) + ky{<7y ® T^) 



(3) 



where rj^] stands for the coordinate at i*'' zigzag chain. It 
can be shown that the valley near fc = 37r/2a in Fig. 1(b) 
originates from the iir_|_-point, the other valley at fc = 
-37r/2a from K.-point^i^i^ 

Since the momentum difference between two valleys is 
rather large, Afc = fc+ — fc_ = 47r/3a, only short-range 
impurities (SRI) with a range smaller than the lattice 
constant causes inter-valley scattering. Long-range im- 
purities (LRI), in contrast, restrict the scattering pro- 
cesses to intra-valley scattering.^ 

In the graphene systems, pairs of time-reversed states 
are formed across the two valleys (Dirac points). In the 
absence of inter- valley scattering for LRI, this ordinary 
TRS becomes irrelevant, while the pseudo time-reversal 
symmetry with respect to the operator T = —i{ay'^TQ)C 
(C: complex conjugation) appears, where the A-B sub- 
lattices act as pseudospin. This corresponds to the time- 
reversal operation restricted to each valley. The bound- 
ary conditions which treat the two sublattices asymmet- 
rically leading to edge states give rise to a single special 
mode in each valley. Considering now one of the two val- 
leys separately, say the one around fc = fc+, we see that 
the pseudo TRS is violated in the sense that we find one 
more left-moving than right-moving mode. Thus, as long 
as disorder promotes only intra-valley scattering, the sys- 
tem has no time-reversal symmetry. On the other hand, 
if disorder yields inter-valley scattering, the pseudo TRS 
disappears but the ordinary TRS is relevant making a 
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FIG. 5: (a) Distribution of the transmission eigenvalues A: 
p(A), at £ = 0.5 for L/a = 50,250,1000, with d/a = 2.0. 
E = 0.5 leads to 3 incident channels. 12,000 samples with 
different impurity configurations are included in the distri- 
bution, (b) Distribution of the dimensionless conductance g: 
x{g), at L/a = 1000 for the same parameter set. 



complete set of pairs of time-reversed modes across the 
two valleys. Thus we expect to see qualitative differences 
in the properties if the range of the impurity potentials 
is changed. 



III. ELECTRONIC TRANSPORT PROPERTIES 
A. One-way excess channel system 

We numerically discuss the electronic transport prop- 
erties of the disordered nanographene ribbons. In gen- 
eral, electron scattering in a quantum wire is described 
by the scattering matrix. "^"^ Through the scattering ma- 
trix S, the amplitudes of the scattered waves O are re- 
lated to the incident waves I, 
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Or 



II 
Ir 



t' 



II 
Ir 



(5) 



Here, r and r' are reflection matrices, t and t' are trans- 
mission matrices, L and R denote the left and right lead 
lines. The Landauer-Biittiker formula^ relates the scat- 
tering matrix to the conductance of the sample. The 



electrical conductance is calculated using the Landauer- 
Biittiker formula, 



GiE) = -Tr(«t) = -giE). 



(6) 



Here the transmission matrix t{E) can be calculated by 
means of the recursive Green function method J^ii^i^i 
For simplicity, throughout this paper, we evaluate elec- 
tronic conductance in the unit of quantum conductance 
{e'^/nh), i.e. dimensionless conductance 5 (i?). We would 
like to mention that recently the edge disorder effect on 
the electronic transport properties of graphene nanorib- 
bons was studied using similar approach, 

In the clean limit, the conductance of the zigzag ribbon 
can be given simply by the number of the conducting 
channel. As can be seen in Fig. [Ijb), there is always 
one excess left-going channel in the right valley (K^) 
within the energy window of < 1. Analogously, there 
is one excess right-going channel in the left valley {K- ) 
within the same energy window. Although the number 
of right-going and left-going channels are balanced as a 
whole system, however if we focus on one of two valleys, 
there is always one excess channel in one direction, i.e. 
chiral mode. 

Now let us inject electrons from left to right-side 
through the sample. When the chemical potential is 
changed from E = 0, the quantization rule of the di- 
mensionless conductance (qk^) in the valley of is 
given as 



9K. 



(7) 



where n = 0, 1, 2, • • • . The quantization rule in the K—- 
valley is 



9K. 



1. 



(8) 



Thus, conductance quantization of the zigzag ribbon in 
the clean limit near E = has the following odd-number 
quantization, i.e. 



9 = 9K. 



9K_ 



2n+l. 



(9) 



Since we have an excess mode in each valley, the scat- 
tering matrix has some peculiar features which can be 
seen when we explicitly write the valley dependence in 
the scattering matrix. By denoting the contribution of 
the right valley (-K"-|-) as and of the left valley {K-) 
as — , the scattering matrix can be rewritten as 



o 



L 

\o-rJ 



Ir 



L 

\i-rJ 



(10) 



Here we should note that the dimension of each column 
vector is not identical. Let us denote the number of the 
right-going channel in the valley Kj^ or the left-going 
channel in the valley K— as ric. For example, Uc = 1 
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FIG. 6: L-dependence of (a) the averaged dimensionless conductance and (b) the variance of the conductance for zigzag ribbons 
with A*' — 10, short-ranged impurity potential {d/a = 0.05, inter-valley scattering), uo = 1.0, and rump. = 0.1. 20000 samples 
with different impurity configurations are in the ensemble average. 
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FIG. 7: Fermi energy dependence of (a) the averaged dimen- 
sionless conductance and (b) the variance of the conductance 
for zigzag ribbons with TV = 10 for various impurity potential 
ranges at L/a — 250. 10000 samples with different impurity 
configurations are in the ensemble average. Here, wo ~ 1.0 



and Uin 



■ 0.1. 



a.t E — Eq in the Figlljb). Thus the dimension of the 
column vectors is given as follows: 



di'm{I~£) = 71c 
dim{I^) — Tie 



1, dim{I^) = He, 



(11) 



r dim{0+) = 
\ dim{0^) = 



1, dim{0^) = n.cj 
dim{0^) = ?ic + 1. 



(12) 



Subsequently, the reflection and transmission matrices 
have the following matrix structures. 
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(13) 



(14) 



(15) 



(16) 



The reflection matrices become non-square when the in- 
tervalley scattering is suppressed, i.e. the off-diagonal 
submatrices are zero. 

When the electrons are injected from the left lead of 
the sample and the intervalley scattering is suppressed, 
a system with an excess channel is realised in the K—- 
valley. The scattering geometry is schematically drawn 

in FiglH Thus, for single valley transport, the r and 

r'__ are ric x {nc + 1) and (ric + 1) x ric matrices, re- 
spectively, and t and t'__ are {ric + 1) x {ric + 1) and 

ric X ric matrices, respectively. Noting the dimensions 
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the system guarantees the following relation: 



FIG. 8: Fermi energy dependence of (a) the averaged dimen- 
sionless conductance and (b) the variance of the conductance 
for zigzag ribbons with TV = 10 for various impurity potential 
ranges at L/a — 2000. 3200 samples with different impurity 
configurations are in the ensemble average. In this regime the 
averaged conductance for short-ranged impurities(d/a ^ 1) 



9k^ 
9'k_ 



9K_ 
9K. 



(19) 



is zero. Here, 



1.0 and Uin 



0.1. 



the conductance g — qkj^ + gK_ (right-moving) and 
g' — g'j^ + g'j^ (left-moving) are equivalent. If the 
probability distribution of {Ti\ is obtained as a function 
L, we can describe the statistical properties of g as well 
as g' . The evolution of the distribution function with in- 
creasing L is described by the DMPK(Dorokhov-Mello- 
Pereyra-Kumar) equation for transmission eigenvaluesi^ 
In the following, the presence of a perfectly conduct- 
ing channel in disordered nanographene ribbons will be 
demonstrated with the help of numerical calculation. Re- 
cently Hirose et. al. pointed out that the Chalker- 
Coddington model which possesses non-square reflection 
matrices with unitary symmetry gives rise to a perfectly 
conducting channel. However, systems with an excess 
channel in one direction has been believed difficult to 
realize. Therefore disordered nanographene ribbons with 
LRI might constitute the first realistic example. It is pos- 
sible to extend the discussion to generic multiple-excess 
channel model, where the m-PCCs appeaii^ Such sys- 
tems can be realized by stacking zigzag nanographene 
ribbonSf2S The electronic transport due to PCC resem- 
bles to the electronic transport due to a chiral mode in 
quantum Hall systemi^i^ However, it should be noted 
that the PCC due to edge states in zigzag ribbons occurs 
even without the magnetic field. 



B. model of impurity potential 



of r and , we find that r__r and r'__r __ 

have a single zero eigenvalue. Combining this with the 
flux conservation relation {S'^ S — SS^ — 1), we ar- 
rive at the conclusion that t t^__ has an eigenvalue 

equal to unity. Note that t' t'l does not have such 

an anomalous eigenvalue. This indicates the presence of 
a perfectly conducting channel (PCC) only in the right- 
moving channels. If the set of eigenvalues for t'__t'^__ 
is expressed as {Ti, T2, - ■ ■ , Tn^}, that for t tL_ is ex- 
pressed as {Ti, Ta, • • • , T„^, 1}, i.e. a PCC. Thus, the di- 
mensionless conductance g for the right-moving channels 
is given as 

= £ r, = i + £r„ (17) 

i=l i=l 

while that for the left-moving channels is 

i=l 

We see that gK_ — g'K_ + 1- Since the overall TRS of 



In our model we assume that the impurities are ran- 
domly distributed with a density riimp, and the potential 
has a Gaussian form of a range d 

V{n)= "^^pf"^^"^) (20) 

ro {random) 

where the strength u is uniformly distributed within the 
range |u| < um- Here um satisfies the normalization 
condition: 

{full space) 

UM exp (-r^d^) /(V3/2) = (21) 

Since the momentum difference between two valleys is 
rather large, Afc — — = 47r/3a, only short-range 
impurities (SRI) with a range smaller than the lattice 
constant causes inter-valley scattering. Long-range im- 
purities (LRI), in contrast, restrict the scattering pro- 
cesses to intra-valley scattering.— 

C. disordered nanographene ribbons 

We focus first on the case of LRI using a potential 
with d/a = 1.5 which is already sufficient to avoid inter- 
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valley scattering. Fig[3] shows (a) the averaged dimen- 
sionless conductance and (b) the corresponding variance 
as a function of L for different incident energies, aver- 
aging over an ensemble of 40000 samples with different 
impurity configurations for ribbons of the width = f 0. 
The variance, which describes the fluctuation of the con- 
ductance, is defined as 



var{g) = {g^) - {gf 



(22) 



The potential strength and impurity density are chosen 
to be uo — 1.0 and n.imp. — 0.1, respectively. As a typ- 
ical localization effect we observe that {g) gradually de- 
creases with growing length L (Figl3|). Surprisingly, the 
ribbons remain highly conductive even at the length of 
L ~ 1500a, i.e. more than 350nm in the real system. Ac- 
tually, {g) converges to (17) = 1, indicating the presence 
of a single perfectly conducting channel. It can be seen 
that (g) {L) has an exponential behavior as 



{g) - 1 ^ exp(-L/e) 



(23) 



with ^ as the localization length. In this paper, the lo- 
calization length is evaluated, by identifying exp(lng) = 
exp(— L/^). Here, g — g ~ 1 (g = g) for the system with 
(without) the perfectly conducting channel. 

Interestingly, the variance for the LRI case shown in 
Figl3l(b) has large values and slowly converges to zero to- 
ward the long wire limit. Also, the variance for higher 
energy modes has double humps structure in the short 
wire regime, which is also indicating the unconventional 
behavior. Such suppression of the fluctuation in the dif- 
fusive regime may be attributed to the level repulsion 
from the transmission eigenvalue of the PCC4i Note the 
variance is zero for E — 0.1, 0.2, 0.3, indicating the exis- 
tence of the perfectly conducting channel. Such peculiar 
features cannot be seen for SRI case, see Figinijb) for 
comparison. 

Fig|4] shows impurity range (d/a) dependence of the 
averaged conductance for various Fermi energy. Clearly, 
it can be seen that the FCC develops in zigzag ribbons 
if the range of impurity gets larger than the lattice con- 
stant. However, the development of PCC gets moderate 
if the incident energy lies at a value close to the change 
between g ^ 2n—l and 5 = 2n-|- 1 for the ribbon without 
disorder. Since at least one of subbands at these energy 
points gives zero group velocity, it can be considered that 
the intra-valley scattering becomes stronger. This fea- 
ture is also for example visible in above calculations as 
the deviation from the limit (g) ^ 1 for = 0.4 where 
the limiting value (g) < 1 (Figl3|). This feature can be 
seen clearly in Figl7]and[51 which will be discussed later. 

We performed a number of tests to confirm the pres- 
ence of this perfectly conducting channel. First of all, 
it exists up to L = 3000a for various ribbon widths up 
to TV = 40 for the potential range {d/a = 1.5). More- 
over the perfectly conducting channel remains for LRI 
with d/a = 2.0,4.0,6.0,8.0, and uq = 1.0, Uimp. = 0.1 
and iV = 10. As the effect is connected with the sub- 
tle feature of an excess mode in the band structure, it is 
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FIG. 9: (a) Magnetic field dependence of the inverse local- 
ization length for various incident energies. The filled 
symbols indicate data sets of systems with SRI being rather 
sensitive to magnetic fields, and empty symbols denote data 
sets for LRI which are almost insensitive to the field. The 
magnetic fiux passing through a single hexagon ring is mea- 
sured in units of the quantum fiux {(f>o = ch/e). (b) Enlarged 
view of (a) of LRI data 5000 samples with different configu- 
rations are included in the averages. 



natural that the result can only be valid for sufhciently 
weak potentials. For potential strengths comparable to 
the energy scale of the band structure, e.g. the energy dif- 
ference between the transverse modes, the result should 
be qualitatively alteredf^ 

As a further test we evaluate the distribution of the 
transmission eigenvalues and dimensionless conductance 
for fixed wire length. In Fig[5l^a), the distribution of the 
eigenvalues A of the Hermite matrix, tt\ is depicted for 
various wire lengths. With growing length L a progres- 
sive separation of the transmission eigenvalues emerges 
with a strong peak close to (localization) and at 1 (per- 
fect conduction channel). The distribution of the con- 
ductance g (trace of the transmission matrix Tr(tt^)), is 
depicted in FiglSKb) for samples in the long- wire limit. 
Obviously, g only distributes above 5 = 1 with a singu- 
larity at 1. 

Turning to the case of SRI the inter-valley scattering 
becomes sizable enough to ensure TRS, such that the 
perfect transport supported by the effective chiral mode 
in a single valley ceases to exist. For a comparison, we 
show the ribbon length dependence of the averaged con- 
ductance in Figini Since SRI causes the inter-valley scat- 
tering for any incident energy, the electrons tend to be 
localized and the averaged conductance decays exponen- 
tially, (g) ~ exp(— L/^), without developing a perfect 
conduction channel. 
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In Fig[71 Fermi energy dependence of (a) the averaged 
dimensionless conductance and (b) the variance of the 
conductance for zigzag ribbons with = 10 for vari- 
ous impurity potential ranges d/a — 0.05, 1.0,2.0,4.0 at 
Lja — 250. 10000 samples with different impurity con- 
figurations are in the ensemble average. With decreasing 
the range of the impurity potential, the averaged conduc- 
tance also gradually decreases and the perfect conduction 
disappears for short-ranged impurities of dja = 0.05. As 
we have already briefly mentioned, the sharp dips can 
be observed at energies where the number of conduct- 
ing channels changes. Also, the variances becomes rela- 
tively large at these energies. The drop of the conduc- 
tance close to E=0 indicates that the two valleys are not 
well-separated in fc-space, because the partial flat bands 
have finite curvature for narrow graphene ribbons such as 
N = 10. Therefore, if the ribbon width widens, the per- 
fectly conduction recovers even close to the zero-energy. 
In FigEl similar figure for Lja — 2000. 

In order to demonstrate that the qualitative difference 
between the two regimes, LRI and SRI, is indeed con- 
nected with TRS, we study the effect of magnetic field 
coupling to the electrons through the Peierls phase. For 
the time reversal symmetric situation resulting from SRI 
scattering the magnetic field removing TRS should have 
a stronger effect than for the case of LRI where TRS 
is broken already at the outset. We use the localiza- 
tion length ^ as an indicator. In Fig[9l the field depen- 
dence of the inverse localization length is shown for var- 
ious incident energies (filled symbols for SRI and empty 
symbols for LRI). Indeed the locahzation length displays 
a stronger field dependence than the LRI. Actually for 
LRI even a so-called anti-localization behavior with in- 
creasing field is visible consistent with recent reports on 
graphene ji^i^i^ Note that for i? < 0.4 only a single chan- 
nel is involved in the conductance such that for LRI no 
localization occurs, i.e. =0. 



IV. GENERAL EDGE STRUCTURES 

As we have seen, zigzag ribbons with long-ranged im- 
purity potentials retain a single FCC. This FCC orig- 
inates for the following two reasons: (i) The spectrum 
contains two valleys (two Dirac iir-|--points) which are 
well enough separated in momentum space as to sup- 
press intervalley scattering due to the long-ranged impu- 
rities, (ii) the spectrum in each valley is chiral by pos- 
sessing a right- and left-moving modes which differ by 
one in number, and so scattered electrons can avoid in 
one channel backscattering. Is the zigzag ribbon the only 
nanographene ribbon showing this effect? Here we would 
like to show that such conditions can be satisfied in the 
graphene nanoribbon with more general edge structures 
except for one case, the armchair edge. In the following 
we will first show that only the armchair edge cannot pro- 
duce the localized edge states. Moreover, the spectrum 
of the armchair ribbon overlays the two Dirac If-points 



at the single momentum fc = and displays therefore 
no separation into two valleys unlike the zigzag ribbon. 
We will then show that general edges can produce the 
necessary conditions to observe a PCC in a disordered 
system. 



A. Continuum approach 

We consider graphene at half-filling in order to ex- 
plore their zero-energy edge states. Here from the tight- 
binding model we derive the stationary Schrodinger equa- 
tion for graphene in momentum space. 





ek 



k 



(24) 



where Cfc = -tJ2texp{ik ■ n) (ti = a(0, 1/v^), T2 
a(-l/2, -1/2V3), and xg = a(l/2, -1/2^3)) and 



(25) 



with ipA{k) and ipsik) are the wavefunctions located 
on the A- and B-sublattice, respectively. The spec- 
trum contains the well-known linear Dirac spectrum 
at two nonequivalent momentum K-^ points, K± ~ 
^(±i,^). In FigHni (a) the lattice structure of 

graphene and the definition of Ti{i = 1,2,3) and coor- 
dinates, (b) corresponding 1st Brillouin Zone, and (c) tt 
band structures are shown. 

In a first discussion concerning edge states we perform 
the following transformation: 



0p(fe) 



72 ^#B(fc)-#A(fc) ' ^ ' 



which yields the Schrodinger equation, 

(e+(7^ - ie^a^)^k = E^k, 



(27) 



where ~ (efc±e^)/2. The structure of this equation is 
identical to that of a BCS problem in particle-hole space 
{ipp, ■0ft,), constituting a Bogolyubov-de Gennes equation. 
The diagonal terms, = Re(efe), formally correspond 
to the kinetic energy for particles and holes, and the off- 
diagonal terms, — ie^ = Im(efc), can be considered as the 
pair-potential of a " superconductor" . Thus we rewrite 
the equation as 



(^feCT^ + Afc(7a;)$fc = E'S'k 



(28) 



In the continuum limit we can approximate — 
— (?i^/2m)V'^ — fi, where fj, is the chemical potential, and 
Afc represents the pair-potential of the superconductor, 
Afc = (e/c — e^) /2i = Im(efc). In Fig[Tll the contour plot 
for (a) (b) are shown. The pair potential shows 



line nodes along the momenta ky = ±\^kx and ky 



0. 
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FIG. 10: (a) The lattice structure ol graphene and the definition ol n {i = 1,2,3) and coordinates. Here ti = a(0, l/\/3), 
T2 = a(-l/2, -1/2^3), and rg = a(l/2, -1/2^3). (b) The corresponding 1st Brillouin Zone. K± = ^(±|, ^), T = (0,0), 
and so on. (c) The band energy structures for vr electrons. 




(a) 



FIG. 11: The contour plot for (a) e+ (b) Afe. The dashed line 
means the negative values. 



Thus close to the Fermi energy we approximate the pair 
potential simply by its angular dependence 



A = |A|cos(36'), 



(29) 



where 9 is the angle of the Fermi vector relative to the 
positive fca;-axis. Interestingly, the "pairing" symmetry 
is odd-parity, an /-wave state. 

From this properties we use now the general rules for 
the presence of zero-energy edge states. -^s^ A related 
discussion can be found in ref^S. Considering a clas- 
sical trajectory, the existence of a zero-energy bound 
state requires that the momentum incident to the edge 
and the momentum of specular scattered outgoing tra- 
jectory lie on the Fermi surface on points which have 
a phase difference of tt for the pair potential. FigfT2l 
shows that this condition is satisfied for the zigzag edge, 
i.e. A{kx,—ky) = —A{kx,ky). On the other hand, 
the phase difference is zero for the armchair edge, i.e. 
A(— fc^:,^:;^) = A{kx,ky). Thus we do not expect zero- 
energy bound states in the latter case. The condition 
that there are no trajectories with a 7r-phase difference is 
only satisfied for the armchair edge. 

This conclusion is consistent with the results of tight 
binding model4ii^ Also, similar conclusion has been 
recently given by the approach of k • p equation by 
Akhmerov.^^ The electronic properties of nanoribbons 




(b) 



zigzag, bearded, cove edge 




FIG. 12: Schematic view of the electron scattering in the sur- 
face region of graphene with the analogy to the Andreev re- 
flection, where the electronic states of graphene were mapped 
to the picture of the /-wave superconductors(see text). Cir- 
cles with shading shows the wavenumber dependence of Afc, 
where shading means negative values, (a) The scattering tra- 
jectories at the armchair edge, which convert the momen- 



tum as kx 



-kx, do not connect gap regions of oppo- 



site sign. However, (b) The scattering trajectories at the 
zigzag, bearded and cove edge, which convert the momen- 
tum as ky ^ —ky, connect gap regions of opposite sign which 
lead to a phase shift ir and zero-energy bound states. 



with clean edge based on the k • p equation can be found 
in refi^i^. 

Let us now turn to the original two-sublattice repre- 
sentation and consider the continuum limit in order to 
analyze possible surface bound states. For zigzag edges 
the X-axis is taken to lie parallel to the edge and the y- 
axis parallel to the normal vector. Along the edge we 



10 



choose the momentum kx- Expanding around the K^- 
point we separate the fast oscillating part by introducing 
Py = kx — 27r/\/3a. Then, we write the wave function as 



Ms)(r) = /A(B)(2/)e^'-Vvfe^ 



(30) 



where fA{B){y) describes the slow y-dependence. Con- 
centrating on the energy close to zero we can also take 
the linear momentum dependence along the y-direction 
so that the effective Hamiltonian has the form 



Hab 
Hba 



fA{y) 
fB{y) 



E 



fA{v) 



(31) 



where 



hsA = h*Ag = —te 



1 



Py 



a (Dk 

'V3l 2 

(32) 

with Dk^ = 2 cos{kxa/2). Now we replace py = —idy and 
formulate the differential equations for fA(B){y)- The 
outer-most lattice sites on the zigzag edge belong to one 
sublattice only, say the B-sublattice. Then we take the 
A-sublattice wave function to vanish /^(y) = 0. For the 
zero-energy i? = we analyze the equation for which 
reads now 



dyfB{y) 



^^'^My) = aB{kx)fB{y) . (33) 



a 1 



This equation has an exponentially decaying bound-state 
solution fB{y) = Ce°"^y whose existence condition is 
ci-B{kx) > 0, resulting in \kx\ > 27r/3. As k^ approaches 
27r/3 the bound state extends deeper into the bulk, be- 
cause its extension is given by 



i{kx) 



a l + cos{kxa/2) 
71 1 - 2 cos{kxa/2) ' 



(34) 



Analogous conditions can be found for other related 

edges, such as the bearded edge where /b(j/) = 0. Then 
we obtain zero-energy bound states for a a = — as > 
which leads to the condition \kx\ < 27r/3 consistent with 
numerical calculations for the tight-binding model. 

Now we turn to the armchair edge which extends along 
the y-axis. Both and K— are projected on the same 
point in fcj,-momentum space. Analogous to the zigzag 
edge we expand the wavefunction around one of the two 
points, say K— = ^ (|, O) and write 

MB)ir) = fA(B){xy'^'e'^\ (35) 

where we extract again the fast oscillating part along 
the X-axis normal to the edge. We linearize again the 

spectrum close to zero energy and take Px = k^ — 47r/3a. 
Then, for fA{B)i the Hamiltonian has the form 



hAB 

hsA 



.fA{x) 
fB{x) 



= E 



fA{x) 

fB{x) 



(36) 



where 

hBA = h*AB = -te-'^^"-/^ _ 1 



(37) 



where we redefined the unit length as a = a/Sg which 
corresponding to the length of the translational vector 
for armchair ribbons. Using Py = —idy and making the 

ansatz 



fA{x) 

fsix) 



(38) 



The boundary condition for the armchair edge require 
{4>A; (/>_b) 7^ as we will see. For = 0, the A has two 
solutions 



a 4 



- sin ^ 

a 2 



and leads to the wavefunction, 



fAix) 
fB{x) 



(39) 



(40) 



the coefficients are determined to satisfy the condition 
that the dimensionless electronic current normal to the 
edge vanishes, i.e. jx = with 



rB{x)fA{x)+fX{x)fB{x). 



(41) 



This leads to the solution 



V2 ll 



■coshCx. (42) 



It is easy to see by inserting this wavefunction into the 
Schrodingcr equation that the only zero-energy states is 
obtained for fcj, = for which f] — Q. Thus no zero-energy 
bound state exists in the case of the armchair edge. 



B. Spectrum and valley structure of other edges 

We now analyze the energy band structure of graphcne 
nanoribbons with general edge structures based on the 
tight-binding model. We will show that the appearance 
of flat bands is a general features (apart from armchair 
ribbons) and yield a two- valley structure with chiral edge 
states similar to the case of zigzag edges. 



1. Bearded Edge and Gove Edge 

We start our discussion with the most simple edge 
shapes having translational symmetry along the zigzag 
axis, called bearded and cove. Although these two edges 
look rather artificial compared with the zigzag edge, they 
are interesting because they show very much the same 
non-bonding edge localization as the standard zigzag 



(a) (b) (c) 




FIG. 13; (a) The structure of a graphite ribbon with two 
bearded edges, (b) The band structure of a semi-infinite 
graphite sheet with a bearded edge, (c) The band structure 
of bearded ribbon for A'' = 10. 



(a) (b) 




FIG. 14: (a) The structure of a graphite ribbon with a zigzag 
and a bearded edge and (b) the band structure for TV = 10. 



edge. A bearded edge is derived from a zigzag edge by 
adding single 7r-electron hopping bonds on each bound- 
ary site as shown in Fig flsy a). This type of edge was 
first studied by Klein. In FigfTSfb). the band structure 
of a semi-infinite graphite sheet with a bearded edge is 
shown. Interestingly, a partial flat band appears in the 
region of < 27r/3, which is the opposite condition in 
the semi-infinite graphite sheet with a zigzag edge, as we 
had anticipated by our continuum approximation above. 

It is interesting to consider a ribbon having one edge 
of zigzag and the other of bearded shape as shown in 
Fig fTW a). Because for this ribbon |A'a — -^^b| — 1, where 
Na{Nb) means the number of sites belonging to the 
A(B)-sublattice, there is a flat band at = all over 
the 1st BZ, as shown in Fig fnf b). The analytic solution 
of this flat band can be easily understood by the combi- 
nation of two edge states for zigzag and bearded edges. 
In the region of |fc| < 27r/3, the electrons are localized at 
the bearded edge, and in the region of |fc| > 2tt/3, the 
electrons are localized at the zigzag edge. At k = ±27r/3, 
the wavefunctions extend over the whole ribbon width. 
It should be noted that this ribbon is insulating at half- 
filling because the flat band has no dispersion and, thus, 
cannot carry currents. Moreover, there is an energy gap 
between the flat band and next subbands. 

Cove edge is a zigzag edge with additional hexagon 
rings attached. A graphene ribbon with two cove edges 
is shown in Fig fTSf a). In FigfTSTb). the band structure of 
a semi-infinite graphite sheet with a cove edge is shown. 
This case also provides a partly flat band in the region 
of |fc| < 27r/3 like the zigzag edge. 

Interestingly both bearded and cove graphene ribbons 
posses two well-separated valleys in momentum space 
and chiral modes in both due to the partially flat bands. 
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(a) (b) (c) 




FIG. 15: (a) The structure of a graphite ribbon with two 
cove edges, (b) The band structure of a semi-infinite graphite 
sheet with a cove edge, (c) The band structure of bearded 
ribbon for = 10. 

Thus in both case PCC are realized as long as the impu- 
rity potential is long-ranged. 

2. General ribbons with \Na — Nb\ ~ 

We extend our analysis to the electronic spectrum of 
nanoribbons for which the ribbon axis is tilted with re- 
spect to the zigzag axis and keep the balance between A- 
and _B-sublattice sites. In Fig lTCT a). we show the defini- 
tions of coordinates and primitive vectors which specifies 
the geometry of the ribbon. For this purpose we intro- 
duce the two vectors, T = (m, n) = mai + na2 and 
W — (0, Z) — la2, where l,m,n are integers. The pure 
zigzag ribbon corresponds to m = —n and the pure arm- 
chair edge is given by m = n. 

FigfTBTb') and (c) show the energy band structures of 
ribbons with the general edge structures of = (0, 20) 
and (b) T = (-4,3) and (c) T = (-6,5) are shown. 
As we expected, the partially flat bands due to localized 
edge modes appear which break the balance between left- 
and right-going modes in the two valleys. Both examples 
are rather close to the zigzag edge so that the two val- 
leys are well separated. In this case PCC can appear. If 
the geometry of the ribbons deviates more strongly from 
the zigzag condition, the valley structure will become less 
favorable for creating a PCC, as the momentum differ- 
ence between valleys shrinks. It is important to note that 
the extended unit cell along these generalized ribbons re- 
duces the valley separation drastically through Brillouin 
zone folding. The length scale is the new effective lattice 
constant ax along the ribbon. Under these circumstances 
the condition for long-ranged impurity potentials is more 
stringent, d being larger than ay and not a. 

V. UNIVERSALITY CLASS 

According to random matrix theory, ordinary disor- 
dered quantum wires are classified into the standard 
universality classes, orthogonal, unitary and symplec- 
tic. The universality classes describe transport prop- 
erties which are independent of the microscopic details 
of disordered wires. These classes can be specified by 
time-reversal and spin rotation symmetry (Table [l|. The 
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FIG. 16: (a) The primitive vectors for nanoribbon with the general edge structures. The translational vector is defined as 
T = (m,n) — mai + na^, and the ribbon width is defined by the vector W — (0,1) — la^- The number of carbon atoms in 
the unit cell is 2{l + l)m. The corresponding energy band structures of W — (0, 20) for (b) T — (—4, 3) and (c) T = (—6, 5). 
Here ot is the effective lattice constant which is given as |T|. 



TABLE I: Universality class 
Universality TRI SRI 
Orthogonal Yes Yes 
Unitary No irrelevant 
Symplectic Yes No 

orthogonal class consists of systems having both time- 
reversal and spin- rotation symmetries, while the uni- 
tary class is characterized by the absence of time-reversal 
symmetry. The systems having time-reversal symmetry 
without spin-rotation symmetry belong to the symplec- 
tic class. These universality classes have been believed 
to inevitably cause the Anderson localization although 
typical behaviors are different from class to class. 

Recently, the presence of one perfectly conducting 
channel has been found in disordered metallic carbon 
nanotubes with LRL— The PCC in this system origi- 
nates from the skew-symmetry of the reflection matrix, 
*r ~ — r,^* which is special to the symplectic symme- 
try with odd number of channels. The electronic trans- 
port properties such system has been studied on the ba- 
sis of the random matrix theoryjS&iJ'^ On the other hand, 
zigzag ribbons without inter-valley scattering are not in 
the symplectic class, since they break TRS in a special 
way. The decisive feature for a perfectly conducting chan- 
nel is the presence of one excess mode in each valley as 
discussed in the previous section. 

In view of this classification we find that the univer- 
sality class of the disordered zigzag ribbon with long- 
ranged impurity potential (no inter-valley scattering) is 
the unitary class (no TRS). On the other hand, for short- 
range impurity potentials with inter- valley scattering the 
disordered ribbon belongs to the orthogonal class (with 
overall TRS). This classification is compatible with the 
magnetic field dependence of the localization length ^ as 



Graphene 

Orthogonal Symplectic 

\ ° 

Nano-graphene W«L^ 
Zigzag 

Orthogonal Unitary 
\ *■ Q 

Armchair 

Ort hogonal ^ Orthogo jal^ 

FIG. 17: Summary concerning for the universality crossover. 
With increasing the range of the impurity potential, graphene 
is known to be the orthogonal for SRIs and the symplectic 
class for LRIs. However, zigzag nanoribbons are unitary class 
for SRIs. Armchair ribbons are irrelevant for the range of the 
impurity. is the phase coherence length. W is the width 
of graphene ribbons. 

shown in FigO Consequently we can observe a crossover 
between two universality classes when we change the im- 
purity range continuously. 

Analogous symmetry considerations can be applied to 
armchair ribbons. In this case the two valleys merge 
into a single one at A; = 0. TRS is conserved irrespec- 
tive of the impurity potential range, if there is no mag- 
netic field. Consequently, disordered armchair ribbons 
belong always to the orthogonal class and do not provide 
a perfectly conducting channel. In view of the fact that 
graphene is known to be symplectic (orthogonal) for LRl 
(SRI))^ it is quite intriguing to realize that the edges in- 
fluence the universality class, as long as the phase coher- 
ence length is larger than the system characteristic size of 
the nanographene system. In the Fig llTl the summary of 
the argument is visualized. In the nano-graphene ribbons 
with general edges, since the two Dirac points are sepa- 
rated in the momentum space with the order of 1 1 /T| , the 
characteristic length causing the crossover is ~ |T|. Thus 
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the PPC is expected to appear even in generic nanorib- 
bons if the sample has only very slowly varying charge 
potential and the intervalley scattering is suppressed. 

VI. CONCLUSION 

The unusual energy dispersion due to their edge states 
gives rise to the unique property of zigzag ribbons. Con- 
cerning transport properties for disordered systems the 
most important consequence is the presence of a perfectly 
conducting channel, i.e. the absence of Anderson local- 
ization which is believed to inevitably occur in the one- 
dimensional electron system. The origin of this effect lies 
in the single- valley transport which is dominated by a chi- 
ral mode. On the other hand, large momentum transfer 
through impurities with short-range potentials involves 
both valleys, destroying this effect and leading to usual 
Anderson localization. The obvious relation of the chiral 
mode with time reversal symmetry leads to the classifi- 
cation into the unitary and orthogonal class depending 
on the range of impurity potential. Since the inter- valley 



scattering is weak in the experiments of graphene, we 
may assume that these conditions may be realized also 
for ribbons. Naturally defects in the ribbon edges and 
vacancies would be rather harmful for the experiment 
making this type of experiment very challengingJ^ii^ 
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